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Abstract. Conventional descriptions of transverse waves in an elastic solid 
are limited by an assumption of infinitesimally small gradients of rotation. 
By assuming a linear response to variations in orientation, we derive an exact 
description of a restricted class of rotational waves in an ideal isotropic elastic 
solid. The result is a nonlinear equation expressed in terms of Dirac bispinors. 
This result provides a simple classical interpretation of relativistic quantum 
mechanical dynamics. We construct a Lagrangian of the form ££ — — S + U + 
K = 0, where S is the total energy, U is the potential energy, and K is the 
kinetic energy. 
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1. Introduction 

The ideal elastic solid has been an important model in the history of physics. It is a 
good approximation for many processes in condensed matter. It was also the basis 
of early theories of light because of its ability to support transverse waves [15] . 
Aether models fell into disuse at the end of the 19 t?1 century when experiments 
failed to detect variations in light speed relative to the direction of the earth's 
motion through space [7]. However, de Broglie's wave hypothesis for matter [5] 
predicts a null result for such aether-drift experiments whether or not the waves 
propagate through a material medium. 

It was recently proposed that a lattice of elastic cells would yield fermions 
and gauge fields equivalent to the Standard Model of particle physics [TT]. This 
mechanical model is similar to the one James Maxwell used to develop the equa- 
tions of electromagnetism [T3J . While a homogeneous elastic solid may be regarded 
as a continuum limit of such a lattice, it is fundamentally simpler in structure. 

It is current practice to describe physical events in terms of discrete elemen- 
tary particles rather than continuous waves. However, the Dirac and Klein-Gordon 
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equations which are used to compute the correlations between particle events also 
determine the wave-like evolution of continuously distributed physical quantities 
such as momentum, energy, and angular momentum densities 12 , 5J. Furthermore, 
the first-order Dirac equation may be regarded as a representation of a second- 
order wave equation with the bispinor wave function encoding the first derivatives 
[H [lOl [9l [3] . Half- integer spin is attributable to the fact that waves propagating in 
opposite directions are independent states separated by 180-degree rotation. Thus 
it is quite possible that a single wave equation may determine the evolution of 
dynamical quantities everywhere in space. 

An ideal homogeneous, isotropic, elastic solid (with no dissipation and con- 
stant elasticity) is the simplest possible model for wave-like processes in three 
spatial dimensions. Hence the physics of an ideal elastic solid should be the foun- 
dation for any analysis of Lorentz-covariant processes in three-dimensional space. 
Yet this simple model has never been properly understood. 

Prior analyses of the elastic solid rely on approximations of infinitesimally 
small rotations or small rates of variation of rotations [6] . Analysis of first-order 
stress and strain yields an elastic wave equation for displacement a(r,£): 



This equation has two obvious flaws. First, since the displacement field is 
defined with respect to positions in space rather than particular elements of the 
solid, the value of the field is transported with the solid as it moves. Inclusion of 
convection by velocity u and rotation by vorticity w yields: 

pd1&. + pu ■ V<9 t a - pw x d t a = (2/z + A) V [V • a] - x [V x a] (2) 

The second flaw is more problematic. A rigid rotation by angle yields a 
non-zero divergence of the displacement field (V • a = 2 (cos — 1)). If varies 
with position either radially or along the rotation axis, then the motion preserves 
volume and the elastic response should be proportional to the shear modulus \i. 
Instead, the term (2/x + A) V [V • a] fails to distinguish between compressional and 
rotational origins of divergence of the displacement field. This is why the equation 
is only valid for infinitesimally small gradients of rotation. 

The goal of this paper is to formulate an exact description of rotational 
(or transverse) waves in an ideal elastic solid. We will accomplish this task by 
describing incompressible motion in terms of rotations rather than displacements. 

2. Basic Assumptions 

We make the following basic assumptions: 

1. The elastic solid is characterized by an inertial density p and coefficient of elas- 
ticity /x, with characteristic wave speed c = yj p/p. 




2. There is a linear response to variations of orientation angle relative to equi- 
librium. This means that an initial static perturbation (with velocity u (r) = 0) 



pd?a = (2/i + A) V [V • a] - /iV x [V x a] 



(1) 
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would yield the response: 

c> t 2 = c 2 V 2 © (3) 

3. The velocity field u has no divergence V • u =0. Therefore the velocity may be 
written as the curl of a vector field: 

u=l[VxJ] (4) 

The vector field J is called the conjugate angular momentum density. It differs 
from the usual definition of angular momentum density J = r x pu in that it 
is independent of the choice of origin and can have arbitrary direction. If u falls 
to zero sufficiently rapidly toward infinity, then integration by parts yields the 
alternative expressions for kinetic energy: 

K = 2 / ^ = 2 / ' ^ 

where w = V x u/2 is the angular velocity, or vorticity. Hence J is the variable 
conjugate to angular velocity for a Lagrangian which depends on u only through 
the (positive) kinetic energy. 

4. We assume that there are no velocity-dependent forces (such as frictional damp- 
ing) . Hence velocity only enters the equations of evolution through the convection 
and rotation of fields. We neglect any effect of the symmetric tensor (piiij + djui) 
on the orientation angle or angular momentum direction. 

Additional assumptions will be introduced in order to simplify the mathe- 
matics, and these may limit the generality of the results. 



3. Equation of Evolution 

Starting from (J3J), we define an angular potential Q such that: V 2 Q = — 4p0. The 
static condition for Q is: 

V 2 {<9 t 2 Q - c 2 V 2 Q} = ( if u (r) = everywhere) (6) 

Define the spin angular momentum as: 

S = d t Q (7) 

When motion is present, it contributes to the time derivative dtS only through 
convection (— u ■ VS) and rotation (w X S): 

V 2 {<9 t S-c 2 V 2 Q + u- VS-wx S} =0 (8) 

From here on, we will consider only wave-like solutions satisfying: 

9 t S-c 2 V 2 Q + u- VS-wx S = (9) 

For oscillatory solutions to this equation, the first two terms are always in phase 
(<9 2 Q — c 2 V 2 Q), whereas the velocity-dependent terms may have different phase. 
However, if the velocity-dependent terms do not add to zero then they must have 
the same temporal phase as the linear terms: 
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u • VS - w x S = n 2 (r) Q (10) 

where ft 2 (r) is some function of position (more generally, fl 2 (r) could have differ- 
ent values for each component of Q. Substitution yields: 

d 2 Q - c 2 V 2 Q + ft 2 (r) Q = (11) 

If Vi 2 (r) is constant and positive, then this is the Klein-Gordon equation. 

Now our only remaining task is to solve for the velocity in terms of other wave 
variables. To do this we note that the above equation ([9]) can be written in terms of 
a four-component complex Dirac bispinor (jp) using the following identifications: 



Sj — dtQj 
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:[V-Q] = -- tytyty] 



c 2 {V x V x Q} j = --ce, 



(12) 



This identification between bispinors and classical variables was made previously 
[3] , where it was also assumed that the velocity could be derived from the spin. In 
this paper we will derive an expression for velocity from the conjugate momentum. 
The matrices tjj are the Dirac spin matrices: 
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(13) 



The matrices cr^dj are the Dirac velocity matrices. One representation for 7 s is: 
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(14) 



The above identifications provide seven constraints on the eight free parameters 
of the Dirac bispinor. In terms of bispinors, the rotational wave equation ([9]) is: 



u • V [^ t c :) '(/'] - eujWk [ip^aiip] = 



(15) 
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Expanding the derivatives yields: 



id 



d t ip + cry a ■ Vip + u • Vip + w • —ip 







(16) 



where (c.c.) represents the complex conjugate. The Hermitian conjugate wave func- 
tion may be regarded as an independent variable (the independent real and imag- 
inary parts of the wave function are linear combinations of elements of ip and tp )■ 
Validity for arbitrary ip* requires the terms in brackets to sum to zero. This yields 
the Dirac equation: 



d t ip + cj 5 u ■ Vip + u • Vip + iw • — tp + i\ip = 



(17) 
(18) 



where x may be any operator with the property: 

Re (ip'cijixip) = 

Since x nas n0 effect on the original equation for Q, we assume it to be zero. 

Now we construct a Lagrange density Jzf '. Lagrange's equations of motion for 
a field variable ip are: 



d[d t iP] 



^ j d[djip] dip 



(19) 



A similar equation holds with ip* replacing ip. A possible Lagrange density for 
rotational waves is therefore: 

en 



if 



iip^d t ijj + ip' -ic7 5 ff • Vl ip + ip> 



-iu • V + w 



iP 



(20) 



Using the Hermitian conjugate of (fT9|) yields simply d^/dip^ — 0. Since the 
complex conjuate of the equation of evolution is also valid, we can require the 
Lagrange density to be real: 

on 



Jzf = Re < - 



i-hp*d t ip + V' t [-ic7 5 cr • V] ip + ip* 



— iu • V + w • 



ip} (21) 



The conjugate momentum to the field ip is p^: 



iip* 



(22) 



d [dtip] 

We assume that we can neglect boundary terms in the integration by parts of: 

J dr 3 w • S = i J dr 3 [V x u] • S = ^ J dr 3 u ■ [V x S] (23) 
The conjugate momentum for r is: 



du 



= Re (-ip^iVip) + -V x ip^tp = pu 



(24) 



This conjugate momentum was derived under the assumption that u is an inde- 
pendent variable. Identification of p r with pu is justified by the lack of external 
forces, implying that pu enters the Lagrangian only through the kinetic energy. 



G 



R A Close 



Making u a function of ip introduces a factor of 1/2: 

if = Re {-i^d t ip + ^ [-ic 7 5 ff • V] V} 
1 

+ M 

Variation of this Lagrange density determines the evolution of rotational waves: 



Re (-V^iW") + -V x ip^-ip 



(25) 



d t tjj + cj 5 g ■ V ip + u • Vi/> + iw • — ip = 



(26) 



where velocity u and vorticity w = V x u/2 are functions of ip as determined from 



4. Dynamical Variables 
4.1. Angular Momentum 

Recall that the velocity is the curl of an origin-independent angular momentum: 

u = — V x J (27) 
2p 

We can separate the conjugate angular momentum (J) into orbital (L) and spin 
(S) components by identifying the orbital (ul) and spin (us) contributions to 
velocity: 



u L = -Re (--0 t iVV') = —V x L 
P 2p 

«s4 vx ^4 vxs ( 2s ) 

Zp 2 Ip 

The incompressibility condition V • u = places an additional restriction on the 
wave function: 

V • pu L = iv • [-ipHVip + i [W>] = X - {-^V 2 ^ + [VV] ^1=0 (29) 

The usual definition of angular momentum J = L + S depends on the choice of 
origin. If we define a rotational velocity as ur = w x r and note that the vorticity 
is the instantaneous angular velocity (w = dt&), then from (19) the conjugate 
angular momentum would have the same form as found in relativistic quantum 
mechanics: 

P ^^0l-^{- rXiV+ 2}^ L + S ( 30 ) 



Exact Description of Rotational Waves in an Elastic Solid 



7 



4.2. Energy and momentum 

In the Lagrange density defined in l|25p. the velocity-dependent term is clearly 
the kinetic energy density. This observation suggests that the Lagrangian has the 
form: 



if = Re j-iVMV* + ^ [~ic7 5 (T • V] 



1 . 



= -£ + U + K 



(31) 



where S = Re (itp^dtip) is the total energy, U = Re (— ip* [ic7 5 cr ■ V] ip) is poten- 
tial energy, and K is kinetic energy density. We will simply adopt these definitions 
without further justification. 

The Hamiltonian is the negative of the energy: 



^ = p^dtip - 5? = Re {ip 



U 



ic7 5 a • 



-ou 



Hamilton's equation for the wave function is: 



d t ip 



cj 5 u ■ 



dp^ d [—iipi] 

We can also define a Hamiltonian operator with dtip 
convention from quantum mechanics): 



(32) 



(33) 



iHip (note opposite sign 



H = ic7 5 ff- V 



1U 



V-w- 



(34) 



The Hamiltonian is a special case (Tg 



of the energy-momentum tensor: 
d^j - (35) 
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The dynamical momentum density P is derived from the orbital part of the angular 
momentum: 

T° = P t = -i^dii, (36) 
This is identical to the dynamical momentum of relativistic quantum mechan- 
ics. Unlike the conjugate momentum, it does not include the rotational motion 
associated with the spin angular momentum. 

The sign of the Hamiltonian (and Lagrangian) is simply a convention. We 
chose signs to make the conjugate momentum parallel (rather than anti-parallel) 
to velocity, thus requiring the Hamiltonian to be the negative of the energy. In 
analogy with plane waves, the function cos (ojt — kx) with lu > is equivalent to 
cos (uit + kx) with oj < 0. We could change the sign of the Lagrangian and Hamil- 
tonian and still preserve the sign of the momentum by using covariant derivatives 
(d/j = (dt, —di)). However, such a relativistic construction is unnecessary and per- 
haps misleading since we are in fact dealing with Galilean space-time (the laws 
of relativity apply to the space of measurements, not to the absolute space-time). 
The energy-momentum tensor components T® — (— S, Pi) may still be regarded 
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as a covariant vector whose magnitude is the scalar S" 1 — P 2 . The equation of 
evolution is of course unaffected by the choice of sign of the Lagrangian. 

5. Discussion 

The description of rotational waves derived here provides a simple classical in- 
terpretation of Dirac bispinors and quantum mechanical operators. Due to the 
nonlinear nature of the equation of evolution, we expect that an elastic solid has 
a discrete spectrum of fermionic soliton solutions. If each soliton is assigned its 
own independent wave function, then cancellation of interference between these 
independent 'particles' would yield the Pauli exclusion principle and require the 
introduction of interaction potentials. Characterization of these interaction poten- 
tials is beyond the scope of this paper. 

Correlations between rotational soliton waves would be computed in the same 
manner as in relativistic quantum mechanics. The reason is that the bispinor wave 
functions, and not the measurement values, satisfy a first order wave equation 
which may be used to directly compute variations in space and time. Therefore 
Bell's Theorem [Tj is not applicable. 

The quantum mechanical Dirac equation is also equivalent to a deterministic 
equation for the evolution of spin angular momentum density. Comparison of the 
classical and quantum equations for spin density may yield new understanding of 
the physical properties of elementary particles. 

6. Conclusions 

A nonlinear Dirac equation describes evolution of rotational waves in an ideal 
isotropic elastic solid. The equation is derived by assuming a linear response to 
variations in orientation and incorporating the influence of motion through convec- 
tion and rotation. The Hamiltonian is represented as a sum of potential and kinetic 
energy. This result provides a simple classical interpretation of Dirac bispinors and 
the associated dynamical operators of relativistic quantum mechanics. 
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